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Abstract 

o, 

^D I We perform the analysis of the existing inclusive deep inelastic scattering (DIS) 

data within NNLO QCD approximation. The parton distributions functions (PDFs) 
2P' and the value of strong coupling constant as(Mz) = 0.1143±0.0013 (exp) are obtained. 

The sensitivity of the PDFs to the uncertainty in the value of the NNLO corrections 
to the splitting functions is analyzed. It is shown that the PDFs errors due to this 
uncertainty is generally less than the experimental uncertainty in PDFs through the 
region of x spanned by the existing DIS data. 
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1. The account of higher-order corrections in an analysis based on the QCD perturbative 
expansions is very important. For the relevant processes measured to the moment the typical 
value of the strong coupling constant Og is 0(0. 1) and the convergence of series in a^ is 
slow. For the deep-inelastic-scattering (DIS) process this problem is especially important 
for the largest and the lowest x regions, where the coefficients of the series contain the 
terms proportional to "large logarithms" . Meanwhile due to great technical difficulties the 
progress in calculation of the higher-order QCD corrections is not so fast. In particular for 
the case of the inclusive DIS structure functions only the two-loop QCD corrections have 
been calculated completely ^. The three-loop (NNLO) case coefficient functions are known 
exactly [^ , while for the corrections to the splitting functions only the even Mellin moments 
up to 8 and some asymptotes were known to the recent time 0, ^. 

An attempt to combine all available information about splitting functions in order to 
obtain reasonable approximation to the exact expressions was done in Refs. P, |[ . The result 
of this study is the set of approximate NNLO sphtting functions in the x space supplied by the 
estimate of their possible variation due to effect of the highest moments. These approximate 
splitting functions have been used in the analysis of Ref. aimed to estimate the effect 
of the NNLO QCD corrections on the shape of the parton distributions functions (PDFs) 
extracted from the global fit. Meanwhile the gluon distribution obtained in this analysis 
turned out to be sensitive to the uncertainties of the NNLO splitting functions given in 
Refs. P, ^. In particular at x ~ 10~^ and Q^ = 20 GeV^ the error on the gluon distribution 
due to this uncertainty is about 35%, which is much larger than the experimental error on 
the gluon distribution obtained in the two- loop analysis of the existing DIS data §]. 

Fortunately the Mellin moments of the splitting functions up to 12 were calculated re- 
cently in Ref. that allowed to elaborate new set of the approximate splitting functions 
with much narrower uncertainty range [|l^. In this paper we describe the results of our 
analysis of the existing DIS data with account of the NNLO QCD corrections. The analysis 
is based on the recent splitting functions given in Ref. 0]. Our main aim is to study the 
effect of NNLO corrections on the PDFs and the value of as extracted from the data with 
a particular attention paid on the errors due to the remaining uncertainty of the NNLO 
splitting functions. 

2. Our theoretical ansatz and the fitting procedure are the same as in our previous 
analysis of Ref. |^, except that now we use the NNLO QCD approximation both for the 
splitting functions and the leading twist (LT) coefficient functions of DIS. Impact of the 
NNLO corrections to the coefficient functions on the values of F2 and Fl is illustrated in 
Fig.|l|, where we give the ratios of these structure functions in the NNLO approximation to 
the ones in the NLO approximation. The NNLO contributions to the coefficient functions in 
the form given in Ref. were used in the calculations. The input for both NLO and NNLO 
calculations was chosen the same as in Ref. |^: The gluon distribution xG{x) = x~'^'^'^{l—x)^, 
the total singlet distribution xS(a;) = 0.6x~°'^(l — a;)'^'^(l + 5x^'^), the number of flavors 
Nf = 4, and the value of Og = 0.2. The largest effect of the NNLO corrections to the 
coefficient functions on the values of F2^l is the rise of Fl at large x. Nevertheless for the 
analysis of existing data this rise is not so important due to sensitivity of the data to variation 
of Fl at large x is rather poor. Much more important is suppression of the structure function 
F2 by ~ 5% at small x since the precision of existing data on F2 is 0(1%) in this region. 
Effect of the NNLO corrections on the splitting functions is demonstrated in Fig.^ where 
the ratios of the logarithmic derivatives of the gluon and the singlet distributions calculated 
in the NNLO and the NLO approximations are plotted. We used in the calculations the 
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Figure 1: The ratios of the leading twist structure functions F2_l calculated in the NNLO 
and the NLO approximations. 
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Figure 2: Ratio of the logarithmic derivatives of the gluon distribution G' = d\nG/d\nQ 
calculated in the NNLO and the NLO approximations (left); the same for the singlet distri- 
bution (right). The dotted curves correspond to the choice A and the dashed curves - to 
the choice B for the splitting functions approximations of Ref. []TU|. 



Table 1: The numbers of data points (NDP) and the x^ values for the separate experimental 
data sets used in the analysis. 



NDP 



Experiment 


proton 


deuterium 


xVNDP 


SLAC-E-49A 


58 


58 


0.64 


SLAC-E-49B 


144 


135 


1.35 


SLAC-E-87 


90 


90 


1.07 


SLAC-E-89A 


66 


59 


1.46 


SLAC-E-89B 


79 


62 


1.12 


SLAC-E-139 


- 


16 


0.57 


SLAC-E-140 


- 


26 


0.90 


BCDMS 


351 


254 


1.17 


NMC 


245 


245 


1.29 


HI (96-97) 


122 


- 


1.12 


ZEUS(96-97) 


161 


- 


1.16 



TOTAL 



1316 



945 



1.13 



approximations of the splitting functions from Ref. |]T0[ and the input distributions from 
Ref. 0. Different curves correspond to the two choices of the splitting functions which give 
the range of the uncertainty of the latter. One can see that the NNLO corrections to the 
splitting functions change the "speed" of evolution moderately: At the scale of Q ~ 10 GeV 
the derivatives change by < 10% at smallest x and even less at the largest x in question (the 
spike at X ~ 0.1 is just due to the QCD evolution has crossover point here and the derivatives 
are very small in this region). As a result the main effect of the NNLO corrections is due to 
corrections to the coefficient function for F2. Figs.|l],0 may be used for the benchmark of our 
NNLO evolution code as well. For this purpose one can compare these figures with Fig. 10 
of Ref. 1^ and Fig. 4 of Ref. [^ correspondingly and convince that the agreement of both 
codes is perfect[} 

The boundary LT PDFs fitted to the data were parameterized within the scheme with 
fixed number of fiavors at Nj = 3. At our starting value of the QCD evolution Qq = 9 GeV^ 
they read 



XUy{x, Qo) = — ^X""(l - X)^"(l + 72X), xdy^X, Qq) 



NJ 



X 



«d/ 



X] 



xus{x,Qo) = ^77ux"-(l-x)''-, xds{x,Qo) = j^x^^^il- 



-XI 



xssix,Qo) = j^ri^x^^^il-xf^' 



xGix,Qo) = Agx'^«(1-x)''g(1+7Pv^+72^x), (1) 

where u, d, s, G are the up, down, strange quarks, and gluons distributions respectively; the 
indices V and 5* correspond to the valence and sea quarks. The parameters A^^, Nj and Aq 
were not fitted, instead they were calculated from the other parameters using the conservation 

"'^The extensive cross-check of the different NNLO evohition codes is underway now and the re- 
sults win be released at the WWW page of the Les Houches workshop "Physics at TeV cohiders" 
( http://pdf.fnal.gov/LesHouches.htin ). 
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Figure 3: The selected PDFs obtained from the NNLO fits with the different choice of the 
NNLO sphtting functions. (Full curves: the la experimental bands for the fit with the 



splitting functions chosen as average of the variants A and B of Ref. [T^; dashed curves: the 
central values for the fit with variant B of Ref. [|l^ ) . 



of the partons momentum and the fermion number. The normalization parameter N^ is also 
calculated from the other parameters in such way that the normalization parameter Ag 
correspond to the total momentum carried by the sea quarks. The parameter ?7s was fixed 
at 0.42 and the other sea distributions parameters were constrained as Ogu = Osd = clss, 

bss = (bsn + &sd)/2. 

The LT structure functions F2^l obtained using these PDFs evolved using the GLAPD 



equations |]rT| were corrected for the target-mass correction and the high-twist contribution 
as well as it was done in our earher analysis of Ref. p|]. The result was substituted to the 
regular expression for the inclusive DIS cross section, which was fitted to the existing data 
varying parameters of PDFs, the value of as, and the high-twist contributions to F2^h- For 
our nominal fit we use the NNLO splitting functions obtained as the average of the variants 
A and B given in Ref. ||10| . 



We used for the analysis the data on DIS of charged leptons off the proton and deuterium 




Figure 4: The la bands for the gluon distributions obtained in the NNLO (full lines) and 
the NLO (dashed lines) analysis at different values of Q. 



targets. The data set coincides in part with the one used in Ref. | 



. The difference from 
and the ZEUS data of 



that analysis is that now we include in the fit the HI data of Ref. 
Ref. [1^ collected in 1996-97 instead of earlier data of these collaborations. Besides, we drop 
the data from the FNAL-E-665 experiment \H\ since they have no impact on the results of 
the analysis due to large experimental errors. It was also checked that inclusion of the high 
Q^ data of Ref. |jl5| does not decrease the experimental errors in the fitted values of PDFs 
and fts and for this reason they were not included in our analysis. The same is valid for 
the ZEUS data of Ref. flSl with Q^ > 300 GeV^ and we also discarded these data points 
from the analyzed data set in order to escape the region where the corrections due to the 
Z-boson contribution should be taken into account. The data points with Q^ < 2.5 GeV^ 
and a; > 0.75 were cut in order to improve the perturbative stability of the results and to 
minimize the effect of nuclear corrections correspondingly. The resulting data set outlined 
in Table |I| spans the region of x = 5 ■ 10-^ ^ 0.75 and Q^ = 2.5 ^ 250 GeV^. 

The statistic and systematic errors in the experimental data were combined in the mini- 
mized x^ using the covariance matrix approach as well as it was done in our earlier analysis. 
The normalization factors for all experiments excluding the old SLAG ones were also in- 
cluded into the covariance matrix, while for the latter the fitted re-normalization factors 
were introduced (see Ref. |p for details). In Table ^ we give the value of x^ obtained after 
the fit and the contributions of each separate experiment to its total value. One can see that 
the value of x^ reduced to the number of data points (NDP) is about unity for the total data 
set. The deviation of the values of x^/NDP off unity for some separate experiments can be 
attributed to the statistical fluctuations in the most cases. This allows one to conclude that 
in good approximation the data can be described by the statistical model with the Gaussian 
probability functions for all errors including the systematic ones and hence the errors in the 
fitted parameters are also Gaussian distributed. The obtained values of the fitted parameters 
with their errors are given in Table ^. 

To examine the sensitivity of our results to the specific choice of the NNLO splitting 



functions we repeated the fit with the choice B for sphtting functions of Ref. |Ty]. The 
results of this fit are compared to the ones of the nominal fit in Table |^ and in FigJ^. One 
can see that the difference of the PDFs values obtained in the fits with different choices of 
the NNLO splitting functions is largest for the gluon distribution at small x and anyway 
does not exceed the experimental uncertainties in the PDFs through the whole region of x 
in question. 

Table 2: The values of the fitted parameters of the leading twist PDFs and the strong 
coupling constant. 

NLO NNLO 









(A±B)/2 


B 


Valence 












du 


0.709 ±0.027 


0.726 ±0.025 


0.731 ±0.025 




bu 


3.911 ±0.051 


4.023 ± 0.049 


4.016 ±0.049 




ll 


1.06 ±0.35 


1.04 ±0.33 


1.02 ±0.32 




ad 


0.706 ± 0.073 


0.762 ± 0.072 


0.792 ± 0.071 




bd 


4.95 ±0.12 


5.15 ±0.13 


5.18 ±0.15 


Glue 












aa 


-0.145 ±0.019 


-0.121 ±0.022 


-0.082 ±0.022 




ba 


8.2 ±1.3 


9.2 ±1.1 


9.9 ±1.0 




1? 


-3.79 ±0.45 


-3.93 ±0.52 


-4.37 ±0.45 




1^ 


7.7 ±1.7 


8.4 ±1.7 


9.4 ±1.5 


Sea 












As 


0.165 ±0.011 


0.1616 ±0.0091 


0.1614 ±0.0082 




O-sd 


-0.1961 ±0.0048 


-0.2088 ± 0.0044 


-0.2068 ± 0.0042 




bsd 


4.7 ±1.3 


5.2 ±1.2 


5.7 ±1.2 




Vu 


1.16 ±0.11 


1.13 ±0.10 


1.09 ±0.10 




bsu 


10.42 ±0.86 


10.72 ±0.84 


10.57 ±0.83 




aJMz) 


0.1171 ±0.0015 


0.1143 ±0.0013 


0.1146 ±0.0012 



We also performed fit to the same data within the NLO approximation in order to check 
the perturbative stability of our analysis. The comparison of the results of this fit with 
the NNLO ones is given in Table |^. One can see that the main difference between the 
NLO and NNLO results is for the value of Og and for the parameters describing the sea 
and the gluon distributions at small x. Nevertheless, as one can see from Fig.^, even in 
this region the shift of the NNLO gluon distribution as compared to the NLO one is of 
the order of magnitude of their experimental errors in the wide range of Q. The same is 
valid for the sea distribution and other distributions are even less sensitive to the inclusion 
of the NNLO corrections. The N^LO QCD corrections as they were estimated in Ref. |TH 



should have smaller effect than the NNLO ones. The reasonable conclusion based on these 
considerations is that the perturbative stability of the obtained NNLO PDFs is better than 
their experimental uncertainties at x > 0.0001. A particular feature of our analysis is that 
our gluon distribution is positive up to Q ~ 1 GeV^ in the region of x > 0.0001, in contrast 
with the gluon distribution obtained in the analysis of Ref. M . 



Our value of 

afNL°(Mz) = 0.1143 ± 0.0013 (exp) (2) 

is by 2cr lower than the value of 



a^^^^^(Mz) = 0.1166 ± 0.0009 (exp) (3) 



obtained in the analysis of the similar data set performed in Ref. |]T^. We should also un- 
derline other differences of those results with ours. Contrary to the results of Ref. |T^ , we 
observe sizeable decrease of the Og value under inclusion of the NNLO corrections (compare 
afLO(Mz) = 0.1171 ±0.0015 (exp) in our analysis and a^^^{Mz) = 0.1155 ±0.0014 (exp) in 
Ref. |T^). In addition, we do not observe the sharp decrease of the error in a^^'"^ as com- 
pared with the error in af ^°. Among the most probable explanations of these discrepancies 
is the difference in treatment of the experimental data. For example one cannot compare the 
experimental errors in as given in Eqn.(^) and in Eqn.(|^) since the latter does not account for 
the systematic errors in data. Besides, the analysis of Ref. ^^ was performed assuming that 



the contribution of the high-twist terms is zero, while we fitted this contribution together 
with other parameters. Evidently, since the high-twist contribution to F2 and the value of 
Qfs are strongly anti-correlated, this may take effect both on the central value and the error 
in as- Nevertheless for the comprehensive clarification of the differences between our results 
and the ones of Ref. [^ a dedicated analysis is needed and we suppose to do it in future as 



well as the comparison with the earlier NNLO fit to the data on the neutrino DIS structure 
function F^^ (see Ref. |T^ and references therein). 

3. In summary, we performed the analysis of the existing inclusive DIS data within 
the NNLO QCD approximation. The PDFs and the value of strong coupling constant 
as(Mz) = 0.1143±0.0013 (exp) are obtained. The sensitivity of the PDFs to the uncertainty 
in the value of the NNLO corrections to the splitting functions is analyzed and it is shown that 
the PDFs errors due to this uncertainty is generally less than the experimental uncertainty 
in PDFs through the region of x spanned by the existing DIS data. The obtained set of 
PDFs may be used to reduce the higher order QCD uncertainty in the predictions of the 
cross sections of the hard scattering processes in the hadron collisions. In particular this 
may be important for reliable estimation of the K-factor for the Higgs boson production (see 



Ref. |T9[ in this connection). 
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